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1 Introduction 



There is a large class of 2D quantum field theories (QFTs) which can be consid- 
ered as perturbed conformal field theories (CFTs). These theories are completely 
defined if one specifies its CFT data and the relevant operator which plays the role 
of perturbation. The CFT data contain explicit information about the ultraviolet 
(UV) asymptotics of the field theory, while the long distance behavior requires care- 
ful analysis. If a perturbed CFT contains only massive particles it is equivalent to a 
relativistic scattering theory and it is completely defined by specifying the S-matrix. 
Contrary to CFT data, the S-matrix data exhibit information about the long dis- 
tance properties of the theory in an explicit way while the UV asymptotics have to 
be derived. 

A link between these two kinds of data would provide a good viewpoint for un- 
derstanding the general structure of 2D QFT. In general this problem does not look 
tractable. Whereas the CFT data can be specified in a relatively simple way the gen- 
eral S-matrix is a rather complicated object even in two dimensions. However, there 
exists a rather important class of 2D QFTs (integrable theories) where the scattering 
theory is factorized and the S-matrix can be described in great detail. These theories 
admit a rather complete description in the UV and IR regimes. 

In this paper we consider the application of the boundary one-point functions 
and boundary scattering theory to the explicit construction of the classical solutions 
describing the boundary vacuum configurations. As an example of an integrable 
QFT we study the simply laced affine Toda theory (ATT), which can be considered 
as a perturbed CFT (a non-affine Toda theory). This CFT possesses an extended 
symmetry generated by a H^-algebra. We consider the boundary conditions that 
preserve this symmetry. This permits us to apply the "reflection amplitude" approach 
for the calculation of the boundary vacuum expectation values (VEVs) of exponential 
fields in ATTs. The reflection amplitudes in a CFT define the linear transformations 
between different exponential fields corresponding to the same primary field of the full 
symmetry algebra of theory. They play a crucial role in the calculation of VEVs of the 
exponential fields in perturbed CFTs |l|] as well as for the analysis of UV asymptotics 
of the observables in these QFTs The boundary VEVs, in particular, contain 

the information about the boundary values of the solutions of classical boundary ATTs 
equations. The information about the long distance behavior of these solutions can be 
extracted from the boundary S'-matrix. The explicit solutions which are constructed 
in this paper provide us the consistency check of CFT and S'-matrix data in integrable 
boundary ATTs. 

The plan of the paper is as follows: in section ^ we recall some basic facts about 
Toda theories and one-point functions in ATTs defined on the whole plane. In sec- 
tion |]we consider integrable boundary ATTs. The boundary one-point functions are 
used to derive the boundary values of the solutions describing vacuum configurations 
in ATTs. We calculate the classical boundary ground state energies and give the 
conjecture for them in the quantum case. In section | we describe the boundary scat- 
tering theory, which is consistent with this conjecture, and construct the boundary 
state. The vacuum solutions describe the semiclassical asymptotics of the one-point 
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functions of the bulk operators in boundary QFTs. This makes it possible to derive 
the long distance asymptotics of vacuum solutions. These asymptotics determine 
completely the explicit form of solutions which we construct in section |^. The solu- 
tions can be written in terms of r-functions, associated with multisoliton solutions of 
ATTs. In the last section we consider the opposite (dual) limit of the quantum cor- 
relation function of the bulk field with a boundary and check that its boundary value 
coincides with that given by boundary VEV. At the end of the section we consider 
the solutions for non-simply laced ATTs which can be obtained by reduction from 
the solutions for simply laced cases. 



2 AfRne and Non-AfRne Toda Theories 

The ATT corresponding to a Lie algebra G of rank r is described by the action: 
where e,, i = 1, t are the simple roots of G and — Cq is a maximal root: 



eo + E 



i=l 



The fields (p in Eq. (|l|) are normalized so that at = 

(V9a(x) ipb{y)) = -Sab log \x - y\ 



(2) 



(3) 



We will consider later the simply laced Lie algebras A, D, E. 

For real b the spectrum of these ATTs consists of r particles which can be asso- 
ciated with the nodes of the Dynkin diagram (or fundamental representations) of G. 



The masses of these particles rrii {i 



, r) are given by: 



m,- = miyi 



m 



— m? 



(4) 



where h is the Coxeter number of G and vf are the eigenvalues of the mass matrix: 



ah 



+ (eo)'^(eo)'. 



(5) 



The exact relation between the parameter m characterizing the spectrum of phys- 
ical particles and the parameter /x in the action (|l|) can be obtained by the Bethe 
ansatz method (see for example @,||^). It can be easily derived from the results of 
Ref. [O and has the form: 



7r/i7(l + 6^) 



Mg)r(^)r(f) 
2r (i) h 



X 



(6) 
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where as usual 7(2;) = r(2;)/r(l — z); x = b"^/ (1 + 6^) and 



/ r \ i/2/i 

HG)=[UnTj (7) 

with rii defined by the equation (||). 

The ATTs can be considered as perturbed CFTs. Without the last term (with root 
Co) the action (|1|) describes NATTs, which are conformal. To describe the generator of 
conformal symmetry we introduce the complex coordinates z = Xi +ix2, z = Xi — 1X2 
and the vector 

g=(6 + l/6)p, P=|Ea>o«, (8) 

where the sum in the definition of the Weyl vector p runs over all positive roots a of 
Lie algebra G. 

The holomorphic stress energy tensor 

T{z) = -\{d.^Y + Q.dl^ (9) 

ensures the local conformal invariance of the NATT. 

Besides conformal invariance NATT possesses an additional symmetry generated 
by two copies of the chiral iy(G)-algebras: W{G)® W{G). The full chiral W{G)- 
algebra contains r holomorphic fields Wj{z) {W2{z) = T{z)) with spins j which follow 
the exponents of the Lie algebra G. The explicit representation of these fields in terms 
of fields dz^Pa can be found in 0. The exponential fields: 

Vaix) = e"-'^(") (10) 

are spinless conformal primary fields with dimensions A (a) = Q -("-Q) ^ ^jj^g fields 
Va{x) are also primary fields with respect to full chiral algebra W{G). It was shown in 
Ref. that fields Va{x) and Vs(a){x) correspond to the same representation of W{G)- 
algebra if their parameters are related by any transformation of the Weyl group W 
of the Lie algebra G acting on the vector a as follows: 

a s{a) =Q + s{a~Q) , seW. (11) 

This means that the fields K(a) for different s G W are the reflection image of each 
other and are related by a linear transformation. The functions that deflne this 
transformation are known as "reflection amplitudes". Reflection amplitudes play 
an important role in the analysis of perturbed CFTs. In particular they permit to 
calculate the vacuum expectation values of exponential flelds in ATTs: 

G{a) = {expa- ip) (12) 
These VEVs were calculated in Ref. [0 and they can be represented in the form: 



G(a) 



^fcr(i^)r(f)x 
2r (i) h 




exp / —[a^e-^'-F{a,t)] (13) 
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with 



F{a,t) = sinh((l + 

^ sinhtsinh(62t) sinh((l + h'^)ht) ^ ' 

here and below subscript a denotes the scalar product of the vector with a positive 
root a, i.e.: 

Oo = a • a; (a — 2Q)a = (a — 2Q) ■ a. (15) 

This expression satisfies many possible perturbative and non-perturbative tests 
for one-point function in ATT. For example it can be easily derived from Eq.(^ and 
from the equation of motion that the bulk vacuum energy E{b) in ATT is expressed 
in terms of function G{a) as: 

niE{h) = h{l+h'^)liG{hei). (16) 

The values of the function G{a) at the special points hci can be calculated explicitly 
and the result coincides with the known expression for the bulk vacuum energy p: 

E{h) = sm{n/h) ^^^^ 

4sin(7rx//i) sin(7r(l — x)/h) 

where x = &^/(l + &^). 

It is easy to see from Eq.(|T^) that function E{b) (and the function G{a) as well) 
is invariant under duality transformation 6— l/6orx— 1 — x. The same is true 
for the bulk scattering ^-matrix in ATT. This remarkable property of duality is very 
important for the analysis of Toda theories. 

In the limit 6 — * the field (p = bip can be described by the classical equations. 
In particular it can be derived from Eq.(p!3D that the VEV of this field (po = (hip) at 
that limit can be expressed in terms of the fundamental weights uji of G as: 

r 

(po = bipo = ^{\ogni-2\ogk{G))uJi (18) 

1=1 

and coincides with a classical vacuum of ATT (|l]). At this limit we have (see Eq.(^) 
that fi = {mk{G))'^ /ATcb"^ + 0(1) and after rescaling and shifting of the field (p the 
action (|l]) can be expressed in terms of the field (p = bip — (po as: 

A = ^/ d^x [i(a^0)2 + m^ELo^.e^"'^] + 0(1). (19) 

The VEVs given by Eq.(|l3D and scattering theory data were used in Ref. to de- 
scribe the asymptotics of the cylindrical solutions of classical Toda equations following 
from the action given by Eg . ([I9|) . In the next sections we use a similar strategy to 
explicitly find the vacuum solutions in integrable boundary ATT. 
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3 Boundary Toda Theories, One-Point Functions 
and Classical Vacuum Solutions 



In this section we consider simply laced boundary Toda theories defined on the half- 
plane H = {x,y;y > 0). The integrability conditions for classical simply laced ATT 
on H were studied in the paper It was shown there that the action of integrable 



boundary ATT can be written as: 



.^^ i=i 



•Abound J d X 

H 

+fiB I dxY.s.e"''-^/^ (20) 



where either all the parameters Si vanish, corresponding to the Neumann boundary 
conditions: 

dy^{x,0) = 0; (21) 

or they are given by Sj = ±1 and the parameter /i^ is related to the parameter /i in 
the bulk (in the classical case) as: 

/i| = fi/nb\ (22) 
For the Lie algebra Ai (sinh-Gordon model) the integrability conditions are much 



less stringent and the parameters sq and si can have arbitrary values ||TT[. Here we 
discuss this problem for other Lie algebras where the choice of integrable conditions 
is rather restrictive. We consider Toda theories with Neumann boundary conditions 
and in the case when all parameters = 1 (with sg = in the non-affine case). 
Really, these two quite different classical theories in the quantum case are described 



by the same theory and are related by a duality transformation {b 1/6) |TB[, 
The cases corresponding to different signs of the parameters Si are more subtle and 
will be considered elsewhere. 

We start from a consideration of the boundary NATTs described by the action ( pOf ) 
without the last term in the bulk action and with parameter sq = in the boundary 
term. The boundary ATTs will be considered as perturbed boundary CFTs. On the 
half-plane with VT-invariant boundary conditions we have only one VT-algebra. In 
this case the currents Wj{z) should be the analytical continuations of the currents 
Wj{z) into the lower half-plane. In particular, they should coincide at the boundary. 
These conditions impose very strong restrictions on the form of the boundary terms 
in the action. It is rather easy to derive from the explicit form of H^-currents 
that the Neumann boundary conditions (^Tj) preserve l^-symmetry. The boundary 



condition (^2]) whose quantum modified version has the form |[T^ 




(23) 
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(with all Si = l,i = l,...r and Sq = 0) describes the dual theory and, hence, also 
preserves PF— symmetry. 

In the IV-invariant boundary NATT we have two kinds of exponential fields. The 
bulk fields Va{x,y) and the boundary fields Ba{x) are defined as: 

Va{x,y) = exp{a ■ ip{x,y)}; Ba{x) = exp{\a ■ (p{x)}. (24) 

These fields are specified by the same quantum numbers as those for the correspond- 
ing fields (0) defined on the whole plane. In particular, they possess the similar 
properties with respect to the action of the Weyl group of G (see section 2), defined 
by Eg. ([Tl|) . Following the standard lines one can introduce and calculate bound- 
ary refiection amplitudes which relate the fields Ba{x) and i?s(a)(a^) (see Ref. for 
details). 

Boundary refiection amplitudes can be used for the calculation of vacuum expec- 
tation values of the boundary exponential fields in ATTs. One-point functions: 

GB(a) = (exp(ia-¥.))^ (25) 
for ADE ATTs were calculated in llTll and have the form: 



with 



2r (i) h 



exp /f[ye-^*-FB(a,t)] (26) 



F ( f]= f(f) V sinh(6a^t) sinh((&(a - 2Q)^ + h{l + b^))t) 
J^B{a,i) /WZ. sinh2tsinh(262t)sinh((l + 62)/,t) ^ 



where for the boundary conditions (|23D the function f{t) is: 



/(t) = 2e* sinh((l + b'^)t) cosh{bH) (2^ 



and for the dual theory which corresponds to Neumann boundary conditions ( piD we 
should make the substitution f{t) — > /^(t): 

/^(t) = 2e"'' sinh((l + b^)t) cosht. (29) 

It is easy to see from the explicit form of Gb{(i) that in the classical limit (6^0 
with bip is fixed) the boundary VEV <^o,b of the field bip for the Neumann boundary 
conditions coincides with the classical vacuum (^o ([TsD in the bulk. For the boundary 
conditions (|23|) this happens only for Lie algebra A^, where both these values vanish. 
For other cases we can derive from Eqs. (p6| - p8D that difference of VEVs of the field 



bif at the boundary and on the whole plane has a form: 



^ , , x-^ /■ , xsinh((/i — 2pQ,)t) sinht , , 

^^'^ fr'n 7 smha;tcosh(l - x)tsmh/it ^ ^ 



a>0 
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Vector G can be expressed in terms of elementary functions at the semiclassical limit 
a; — 0, at the self-dual point x = 1/2 and in the dual limit x — *■ 1. We discuss the 
dual limit in the last section. Here we note that at a; — > vector has a limit: 



oo 



^ = - E « / - ""^^^\:'f"^'^ anht. (31) 
^ ' t sinh(/it) ^ ' 



Q>0 







For D and E series, which we consider below the numbers 2la = h — 2pa are always 
even. For the following it is convenient to define: 

E,, = exp(i^9-e,)=n.>o(P/J-^'-" (32) 
where pi = 1/p-i can be written in terms of trigonometric functions as: 

, ,°f dtsmh2lt , , -A- cos(7r(/ + 1 - 2A;)//i) , , 

p2 = exp / — — tanht = TT — ^t-ttV^. 33 

t sinhht ' ^\ cos{n{l -2k) h) ^ ' 







In particular: 



E, = hz^EM2zl (34) 

Ph/2-l 

here q{G) = max, n,; q{D) = 2, q{E^) = 3, q{E-j) = 4, q{E^) = 6. 

The vector d is simply related to the boundary soliton solution that describes the 
classical vacuum configuration. The classical problem for this solution (f){y),y > can 
be formulated in the following way. We are looking for a solution to the classical Toda 
equations, which decreases at y — oo and satisfies at y = the boundary conditions 
that follow from the action (|20|) . After rescaling and shifting (see section 2) the field 
= b(f— (fo satisfies the equation: 

r 

dy(f) = m^Yl ^i^i exp(ei -0); no = 1 (35) 

i=0 

and the boundary condition at y = 0: 

r 

dycj) = m^^y/nleiexp^lei ■ 4)). (36) 

i=0 

The vector is equal to the boundary value of this solution and hence, it completely 
fixes the boundary soliton. At ?/ = we have 



0(0) = ^9; 9^0(0) = ^^iEi. (37) 

i=0 

This gives us the possibility to study Eq.(^) numerically. The numerical analysis 
of this equation shows that for these and only for these boundary values a smooth 
solution decreasing at infinity exist (see appendix C). 

It is natural to expect that solution 0(?/) to the Eqs. (P^), (^) can be expressed 
in terms of tau-functions associated with multi-soliton solutions of the classical ATT 
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equations (see, for example |16|,118 and references there). For the D and E series 
of algebras (besides the cases and |]T6| ) the explicit form of these solutions 
was not known. We postpone the discussion of these solutions to section 5. Here we 
consider the classical boundary ground state energy which can be defined as: 



1 y 

= ^[2^E v^^^ + j dy{\{dy^? + EU n.{e^^-^ - 1))]. (38) 

We note that numerical values of defined by Eq.(|31|) are rather small for all 
G, and the functional ( ^8|) can be calculated to a good accuracy using the bilinear 
approximation: 

4tL = £p E v^^^(2 - logi^.) + o{m- 

More careful analysis of Eqs.(|35|),(|36|) (see section 5) gives us the reasons to write 
the expression for boundary ground state energy in the following form. Namely, we 
denote as Sm(G') the sum of the masses of all particles in the ATT: 

T.m{G) = j2mi = mi2^i = m-tr{M^/^). (39) 

i=l i=l 

This has the following values for simply laced Lie algebras: 

^ / . N / , ,N ^ / N 2m cos(7r/4 — 7r/2/i) 
A„ = 2mcot 7r/2/i ; S^p^ = ■ ,nu^ ' 

1/2 cos(7r/8) _ 2msin(27r/9) _ 



T.^{E,) = m(6-2V3)V2 7 Z '^ ^rn{E^ 



sin(7r/2/i)' ""^^^^ sin(7r/2/i) ' 

E„(E8) = 4m(V3 sin(7r/5) sin(7r/30))^/^ w^Vom • (40) 

sm [rr/ 2nj 



The classical boundary ground state energy ( P8| ) can be written in terms of these 
values as: 

4oL{G) = ^ tan(7r/2/.)S„(G). (41) 

In the quantum case the boundary ground state energy Sl'^aund "^i^^ have contri- 
butions coming from the boundary term in the Hamiltonian and from the bulk fluc- 
tuations around the background solution. The contributions of the first type can 
be calculated using the explicit expression for vacuum expectation values (|^). For 
small h the first quantum correction of the second type can be expressed through 
the boundary 5"— matrix at 6 = (see, for example |jT9[,[0). We will discuss these 



boundary S'-matrices below. Here we note that in the strong coupling region 6^1 
our theory is described by the weakly coupled dual ATT with Neumann boundary 
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conditions (|2l|) . In the strong coupling limit the dual theory is a set of r free bosonic 
theories with masses rrii. The boundary ground state energy for the free massive 
bosonic theory with Neumann boundary conditions and mass rrii can be easily calcu- 
lated and is equal to rrii/ 8. The first perturbative correction in the weakly coupled 
dual theory can be also evaluated with the result: 

^ItnAG) = (l + ^ cot{rr/2h) + 0(1/6^)) . (42) 



Both asymptotics 6 — > (^) and b ^ oo (^) are in agreement with the following 
conjecture for boundary ground state energy: 



cii) _ sin(7r/2/i)S^(G') 
boundy ) 8sin(7ra;/2/i)cos(7r(l-a;)/2/i) ^ ^ 

where x = IP' / (1 + IP). 

The nonperturbative check of this conjecture can be made using the boundary 



Thermodynamic Bethe Ansatz equations |21]. The kernels in these nonlinear integral 



equations depend on the bulk and boundary S'-matrices. In particular boundary 
ground state energy can be expressed in terms of the mass rrij of the particle j 
multiplied by the ratio of Fourier transforms of logarithms of boundary ^-matrix Rj 
and bulk amplitude 5*-,^ taken dX uj = i. To check our conjecture and to construct the 
explicit solution of Eq.(B^) we need the information about boundary 5— matrix. 



4 Boundary S-matrix and Boundary State 

The boundary ^-matrix (reflection coefficient) in the ATT for the particle j corre- 
sponding to the fundamental representation vrj(G) can be defined as: 

\3,-e)o^t = Rjm3.e)^n (44) 

where is the rapidity of particle j. 

The refiection coefficients for the ATT with boundary conditions ( ^Sj) were 
conjectured in [P2],p3[. They can be written in terms of function: 



sm ( - + 



in the following way: 



(45) 



R-{e) = Y[{a - l){a - h){-a + x){-a + h + 1 - x). (46) 

a=l 

Unfortunately, we were not able to find in the literature any conjecture for other Lie 
algebras consistent with duality properties discussed above. So we will give here a 
conjecture which naturally generalizes the refiection coefficients (|4^) to other simply 
laced Lie algebras. To do this we rewrite Eq. (EBf) in the form: 
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R,{e)=eM-^m) (47) 

where 

CiCi 

/dt 2h9t 
— sinh( )[sinh((l -x)t)sinh((/i + x)t)A,(A.,t) -2] (48) 
t vr 



and 

. , 8 sinh ?t sinh(/i — ?)t 

AJAr,t) = ^ ; , ^ ■ 49 

The natural generahzation of these equations can be written as: 

R.(e) = $j(^) exp{-i6j{e)) (50) 
where ^j{9) are the CDD factors, satisfying the conditions: 

^j{9)^j{-e) = 1; $(e)$(5 + ivr) = 1 (51) 

and the function 6j {0) is defined by Eq. (|48|) , where we should do the substitution 

/^j{Ar,t) Aj{G,t) with: 

^jiG.t) = -4rr. k(2cosht-2)5^„ + e^-e„)-'l ... (52) 



cosh ht 



n 



The most important part of this conjecture is that for the particles j corresponding 
to fundamental representations TTj{G) (or the node j of Dynkin diagram) with nj in 
Eq.(^ equal to 1, the CDD factors ^j{9) in Eq.(|50D are equal to one. 

This statement fixes completely the boundary S-matrix for Lie algebras Dr,EQ, E-j. 
We denote by Rf{G,6) the fundamental reflection coefficients. This means that 
all other amplitudes Rj{0) can be obtained from the fundamental reflection factors 
Rf{G,9) by application of the boundary bootstrap fusion procedure |TI[|,[Q. It is 
easy to see that Rf{Dr, 6) = Rr{9) = Rr-i{9), where particles r and r — 1 correspond 
to the spinor representations of Dr, Rf{EQ,9) = Ri {6) = Rj{9), where particles 1 
and 1 form the doublet of lightest particles in the Eq ATT and Rf{EY,9) = Ri{9), 
where 1 is the lightest particle in the E-j ATT. For all these three cases the CDD 
factors $/(6') = 1 and the reflection coefficients are given by Eqs.(0),(^) with j = f 
and 

_ 8Binh(M/2).inh((V2 + ,-l)t) 

' ^ sinh(gt) sinh(2/it) ^ ^ 

where q(G) = maxj n^; q{D) = 2, q{EQ) = 3, q{Ej) = 4. 

As an example of the application of the boundary bootstrap equations we give here 
the reflection coefficients Rj for the particles j = 1, 2, n — 2 in the ATT, which 
can be obtained from the amplitude Rj (9) = Rr{9) = Rr-i{9). These functions Rj 
have the form (|50|), where: 
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16 sinh jt cosh(( V2 - j)t) sinh(/it/2) 
^^■(^-') = sinhtsinh2/.t ^'^^ 



and the CDD factor can be written as: 



oo 



$j = exp y J Y sinh(2/i^t/7r) sinh(l - x)tcosh.{xt)iljj{Dr, t)j (55) 

where iljj{Dr,t) = 8sinh(j — l)t sinhjt/(sinh2t cosh/it) 

The Lie algebra Eg has no fundamental representations with rij = 1. The lightest 
particle in this case is associated with the adjoint representation and Rf{Es,6) = 
Radid)- The adjoint representations for Lie algebras D^E have Had = 2 and, hence, 
CDD factors should appear. The reflection coefficients Rad{0) for Lie algebras D and 
E can be written in the form (0), with: 

8cosh((g-l)t)cosh((V2-g)t) 
^"'^^^ ^> - cosh(/.t/2)cosh/it 

where q{G) is defined above and the CDD factor $ad(6') is defined by Eq. (|55|) with 
function 

, ^ , 8 sinh t , , 

The amplitudes Rj{0) following from the scattering theory described above have 
no poles corresponding to the bound states of particles j with boundary. The only 
poles with positive residues which appear in some of the functions Rj{0) are the poles 

at 6* = i7i/2 

^^^^^ - ^^^^ 

corresponding to the particle boundary coupling with zero binding energy. For exam- 
ple, the lightest particle whose boundary amplitude possess such pole is the particle 
vTiad corresponding to adjoint representation. It is easy derive from Eqs. (p6| , |57D that 
corresponding residue has a form: 

DUh) = 2 cot tan tan %T,_,%/,_, (59) 

where % = tan(f - ^) tan(f - ^i^) cot (f - f) cot (f - 

The amplitudes Rj{9) and their residues can be used for the construction of the 
boundary state in integrable QFT. It was shown in Ref. |[Tl|] that this state can be 
written in terms of numbers Dj and functions Kj{6) = Rj{iTi/2 — 9) as: 



\B) = exp 



|0) (60) 
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where Aj{6) are the operators creating the asymptotic states with rapidity 6 i.e. 

Me)\o) = \j,e),^. 

We will be mostly interested by the linear in the operators Aj part in the expansion 
of boundary state \B) 

i5) = [i+E^.m(o)]io) + ... 

which determines "one particle" contributions to the correlation function: 

J',iy) = {0\b^iy)\B)-{0\b^\0) (61) 

where the variable y > plays the role of Euclidean time. These "one particle" 
contributions can be written as: 

J^Uy) = j2D,{b){0\b^mj,0)^nexp{-m,y) + ... (62) 
i=i 

The boundary soliton (j){y) (see section 3) coincides with the semiclassical limit of this 
correlation function i.e. (f){y) = linifc^^o-^b iu)- In the weak coupling limit one particle 
matrix elements can be easily calculated and have a form: 

(0|^(0)|j,0),„ = v^ej + 0(62) 
where C,j are the eigenvectors of mass matrix: M^j = v'j^j-, satisfying the conditions: 

^^■^, = 5i,■, p-e, >o. (63) 

The corresponding "one particle" contributions to the boundary solution 0(?/) will 

be: 

r 

(t>{y) = H dj^j exp(-mj?/) + ... (64) 

where 

dj = \im^bDj{b). (65) 

For example the coefficient dad corresponding to the lightest particle mad in expansion 
(|6^) and defining the main term of the asymptotics aX y ^ oo can be extracted from 
Eq.(pUD and has a form: 

dad = 2^/1 tan ToTq^iTh/2^g (66) 



where = limfe^o % '■ 



^ ,71 7r(fc + l), f 7T 7Tk\ , , 
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The coefficients dj as well as the vector d (or numbers Ej) fix completely the 



solution (j){y) of Eqs.( p5|j36| ). They determine the contribution of the zero modes of 
the linearized Eq.(^) and make it possible to develop in a standard way the regular 
expansion at large distances. If our scattering theory is consistent with conformal 
perturbation theory this expansion should converge to the boundary values at ?/ = 
At the end of this section we note that the analysis of the boundary Thermody- 
namic Bethe Ansatz equations with kernels depending on the refiection coefficients 
written above gives an exact agreement with Eq.(^3]) for the boundary ground state 
energy. 

5 Boundary Solutions 

It is natural to assume that boundary vacuum solutions can be expressed in terms 
of r-functions associated with multisoliton solutions of ATTs equations. This means 
that field can be written as: 

r 

(P{y) = -J2^i^^STi{y); Ti{y)^0,y^Qo (68) 
where functions Ti{y) satisfy the equations: 

- r:\y) + {r^r = (Uy'^ - r^y. (69) 

here Xij is the incidence matrix of the extended Dynkin diagram of G. 

The classical boundary state energy can be expressed in terms of numbers Ei and 
the boundary values of r-functions [jl6[ as: 



Functions Tj corresponding to multisoliton solutions to the ATT equations are given 
by finite order polynomials in the variables 

Zj = exp^—rrijy) = exp^—mujy). (71) 
The general ansatz for these functions has a form: 

n{y) = E - E Ytl,XtiZ,f'-{trZ.)'' (72) 

fci=0 A,v=0 

where Yq^]q = 1 and all other coefficients can be derived from Eq.(|69D if we fix the 
normalization for the linear in variables Zj terms in the expansion (|7^): 

r.(y) = 1 + E Yo%..ohZ, + 0{Z'). (73) 
i=i 

These coefficients should satisfy the condition: 
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-Ee.n%..o = A^^-- (74) 

i=0 

After we fix the factors Nj all coefficients Y^^^ in Eq.(^) are universal numbers 
depending only on G. The parameters tj in the Eq. ( [72D are simply related with 
parameters dj defined by the expansion (|64D. It is easy to see from Eqs. (|6^, [73| , [7^ ) 
and Eq.(|6^) that 

tj = dj/Nj (75) 

The parameters dj can be extracted from the boundary scattering theory described 
in the previous section. This gives us the possibility to derive the explicit solution to 
Eqs.( 



5.1 Dr boundary solution 

The boundary values of this solution can be described by the numbers Ek (see section 
3) which in this case are: 

Eq = El = Ej—l = Er ^ 



h sin (vr/Zi) 



^ ^ cos (nk/h) cos (7r(r -k)/h) plpl_k /^^x 

cos(7r(2A;-l)/2/i)cos(7r(2r-2A;-l)/2/i)' ^ ' 

These numbers possess all the symmetries of extended Dynkin diagram for Dr. The 
same is true for the boundary solution 0(?/). It means that in Eq. ( |68D tq = ti and 
Tfc = Tr-k- For the Lie algebra Dj. the following relation holds: 

r 

— E rij COS (j7r(2i — l)/h)ei = 2\/hsm{nj /h)^j (77) 

i=0 

so that it is natural to choose Nj{D) = 2^/hsm(-Kj /h). It is convenient also to 
introduce the functions: 

To = f 1 = r,,_i = fr = Tg ; ffc = Tfe; k = 2, ...r - 2. (78) 

The particles which are invariant with respect to symmetries of Dynkin diagram 
are particles 21; 21 < r — 2 with masses 17121 = my/S sm{27rl / h) . The analysis of the 
boundary S -matrix (|5l) gives that only amplitudes R2i{0) have a pole at 6* = i7r/2. 
The parameters t2i = d2i/N2i{D) can be derived from the Eqs. (|5^ , |55| ) and have a 
form: 

_ tan (7r/4 - irl/h) ' tan2(7r(2z - l)/2/i) 

~ 2cos2(7r///i) W tan2(7rz//i) ^ ' 

Parameters t2i and functions Z21 = exp(— ymv^ sin(27r///i)) are defined for 21 < r — 2. 
To write the solution in the most short form it is convenient, however, to continue 
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these values to / < r — 2. To continue the parameters to / < r — 2 we can use 
Eg. ([79D ■ In this way we obtain: 



th-2i — —hi', Zh-2i — Z21 (80) 
We introduce also lu = exp(27ri//i). Then the boundary solution can be written as: 



4>{y) = -|ELo^ieilogri(?/) 



^1) 



where: 



1 1 r-2 

(Ti=0 (Tr-2=0fc = l 

^ /sin (vrfm — 



2k 



2(7 m 



2) 



sin (Tc{m + n)/h) J 

As an example of application of the last equation we give here the result for D7 
solution. This solution is characterized by three function tq, 72,73 ([78|) which are: 

\ 2 



T"0,3 



\±^t^Z, + -hZ2±—tan'' 
2 2 4 



TT 
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tlZit2^2 



1"2 



2^2-^2 + 



3(tiZi)2 , (t2Z: 



+ 



TT 



— tanM— (tiZit2^2 
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4 4 16 

We note that in the usual basis of roots of Dr. eo0 = — 0i — (p2', = (pk — 4>k+i, k = 
l...r — 1; = 0r-i + 4>r the solution can be written as 0^ = log (rfc_i/rfc^. It can be 
checked that boundary value of (j){y) defined by Eqs. (|81| , |82D coincides with {}, defined 
by Eq. (|76D and the classical boundary ground state energy (^) can be derived using 
Eq.m. 



5.2 Eq boundary solution 

The boundary values for this solution are: 

E2 = E^ = E^={E^E^Y^'^ (83) 

Here and later we shall use the numeration of nodes j of Dynkin diagram correspond- 
ing to the particles nij (see Fig.s [1-3]) as follows: 



)m2 

O O 6 O O (Fig.l) 

mi ma m4 m-^ mj 
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The solution possesses Z3 symmetry of this diagram. It means that: Tq = Ti = r_ 
and T2 = Ts = r_. The particles that are invariant with respect this symmetry are 
^2 = ""^ad and m^. The analysis of boundary S'-matrix shows that only amplitude 
R2 = Rad{d) possesses the pole aX 9 = m/2 with the residue (pW). In this case only 



one particle contributes to the solution. It is convenient to take N{Eq) = 2\/3h. 
Then the parameter ^2 = dad/N is: 

h = tan ToTsTa/v^ = tan tan tan (84) 
The variable Z2 = exp(— mz/2y), where z/| = 2(3 — ^/S) and the solution has a form: 

To = 1 + (2 + V3)t2^2; ra = 1 - 2t2^2 + ((2 + V3)t2^2) ' ; 

r4 = l-3V3t2^2-3V3(2 + V3)(t2^2)' + ((2 + V3)t2^2)'. (85) 

It is easy to check that 0(0) = and the classical boundary ground state energy ( [70D 
coincides with (ETI). 



5.3 £"7 boundary solution 

The boundary values (|32| ) for this solution are: 



Ps 9 V9/ P3P7 
E, = E, = E, = ^-^- E, = ^^. (86) 

PlP2P3Pm PiPe P1P2P5P7 

Numbers E^ possess Z2 symmetry of the extended Dynkin diagram 



o- 




-o- 



-o 



-o- 



-o o (Fig.2) 



m2 ms mj 1714 mi 



This means that in Eg . (|68D tq = ti, T2 = T4 and Ts = tq and only Z2 even 
particles m2 = mad, ""^4, and my can contribute to this solution. The analysis 
of the boundary S'-matrix (which can be obtained by boundary bootstrap fusion 
procedure from the amplitude Rf{9) (p3D) shows that only amplitudes R2 = Rad, Ra 
and Rj possess the poles aX 9 = in /2. These amplitudes have a form (^) where the 
phases 5j and CDD factors <I>j are defined by the functions Aj(ii^7,t) and il)j{E'j^t) in 
Eqs. (|48|j5^j55| ). For amplitude R2 = Rad these functions are given in section 4. For 
amplitudes R4 and Ry they are: 

16 sinh 5t cosh 4t 
cosh ht ' 
8 sinh 3t sinh 5t sinh 8t 
sinh t sinh 2t cosh ht 



^ 16 cosh t cosh St cosh 4t_ ^ 

cosh(/ii(:/2) cosh /it ' 
^ 8 cosh t cosh 2t sinh 6t ^ ^ 

sinht cosh(/ii(:/2) cosh /it' 
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The coefficients dj can be extracted from these amphtudes using Eq.(|65D. It is con- 
venient to take N{Ej) = 2\/3h. Then the parameters tj = dj/N can be written in 
terms of numbers (^) as: 

V3t2 = tan (^1^) ToTgTs; V^U = ta.n (^^^ ToTiTsUTr, 

V3t7 = tan (^■^^ToT^T2T^T^T,\Te)-\ (87) 
The corresponding variables Zj = exp(— mz/j?/), where: 

z/2 = v/32sin(7r/9)sin(27r/9); 1/4 = y32sin(7r/9)sin(47r/9); 

z/7 = y32sin(27r/9)sin(47r/9). (88) 

We give here the exphcit form for function tq. All other r-functions can be derived 
from it using Eg . (|69|) . For example it holds: 

T2 = - TotI; + (rl^y (89) 

with similar relations for higher r-functions. As an example of application of this 
equation we give the explicit expression for T2 in the form (^) in Appendix A. 

For all simply-laced algebras the number uq is equal to one and function tq in 
Eq.([72D can be written as: 

^0 = E E E Y^^i^^it^z^nuz.ntjz.r m 

<JI=0 (J2=0 O"3=0 

It is convenient to simplify the notations. The variables cTj take the values or 1. 
We can introduce the variable 1 < i < 3 indicating which are equal to one. The 
corresponding coefficients Y^^^^^^,^ we denote as yi,yij,yi23- For example yi = ^100, 
?/i3 = Y^qI and so on. With this notation we have: 

2cos2(27r/9) _2cos2(7r/9) _ 2 cos2(47r/9) 

~ cos(47r/9) ' ~ cos(27r/9) ' ~ cos(7r/9) ' 

8cos^(27r/9) 8cos4(47r/9) ,^ ,^ ,^,,7 

yi2 = r7^'^i3 = ' ' \ yi3= 2cos 47r/9 ^ 

cos(7r/9) cos(27r/9) 

yi23 = (2cos2(47r/9)/cos(7r/9))^ (91) 

It can be checked that 0(0) = "d, where is defined by Eq.(p6|) and the classical 
boundary ground state energy (|70D coincides with that given by Eq.(^). 

5.4 boundary solution 

The boundary values for this solutions follow from Eq.(^), where Eq is given by 
and 

PsPePio PsPioPu PiPaP^Pu 

P1P5P13 PiPsPbPrPu P3PePi2 

p P2P3P5P12P13 p P5P6P12 p P5P9 

-c/4 = ; J^5 = ; -£^6 = ; 



PePiPsPu P1P2P9P11 P2P4PS 

P1P2P7P8P9P13 pmoPnPu 

Ey = ; -t/8 = • l92j 

P4P5P10P10P11 P5P9P12P13 
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The extended Dynkin diagram 

1714 



o o o o o o o o (Fig.3) 

m2 rriQ rrig mj m^, mi Q 

in this case has no symmetries and all particles can contribute to the solution; however 
the analysis of the boundary 5- matrix shows that only amplitudes Ri = Rad, R2, R5 
and Rs have the poles at 6' = i'7i/2. The amplitude Ri = Rad is given in section 4. 
Three other amplitudes with these poles are defined by the functions Aj{EQ,t) and 
tpjlE^jt), which have a form: 

^ 16 cosh 3t cosh 5t cosh 6t _ ^ 16 sinh 7t cosh 6t _ 

cosh(/it/2) cosh /it ' cosh /it ' 

^ 8 cosh 4t cosh 5t sinh 6t ^ 8sinh3t /sinh lit sinh 12t \ 

sinht cosh(/;,t/2) cosh/it' coshht \ sinhtsinh4t y' 

^ 2 sinh 6t sinh lot _ ^ 8 sinh 5t sinh 9t sinh 14t 

sinh^ t cosh ( /it /2 ) cosh ht ' sin t sinh 2t cosh ht 

These amplitudes determine the coefficients dj (see Eq.(^)). It is convenient to take 
N[Es) = Then the parameters tj = dj/N can be written in terms of numbers 

Tk ® as: 

^1 = tan (^^^ ToTsTg; ta = tan (^^j ror3T5T6T9Tii(T8)~^ 
t5 = tan (^-^^ tan TrTniT^TfiTsTioTuy^ 

ts = tan" (^-^^TinT.TjTgTuTMTsT^.T^^yK (93) 
The corresponding variables Zj = exp(— mz/j-y), where: 



1,5 



3(5 - v^)/2 T V3(5 + V5)/2; 



2,8 



3(5 + v^)/2 T V3(25 + liy5)/2. (94) 



We give here the explicit form for function tq. All other r-functions can be derived 
from it using Eg . (|69|) . For example function ti is related with tq by Eq.(|89|). Function 
To can be written in the form: 

O"l=0 (T4 = 

The coefficients Y^'^^^ for this function are given in Appendix B. It was checked with 
the accuracy allowed by Mathematical^ that 0(0) = ^ and the classical boundary state 
energy (ffOD coincide with that given by Eg . (p|) . 

§ Copyright ©Wolfram Research Inc. 
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6 Concluding Remarks 



1. In the previous sections we studied the semiclassical asymptotics of integrable 
boundary ATTs. In particular, we derived the solution (f){y) describing the limit 
6 — > of one-point function J^b{y) of the bulk operator with a boundary. Here we 
shortly consider the opposite (dual) limit 6 — »• cxo of this function. The boundary 
value of function Thijj) is given by Eq.(^), and at this limit is equal 

For example for the Lie algebra this equation can be written as: 



ei • = Cr-i ■ = er ■ = -7T/{2h) cot (27r//i) 
-7r/(2/i)[cot {2'Kk/h) + cot (27r(r - k)/h)]. 



(97) 



At the limit b ^ oo our theory is described by the weakly coupled dual theory, 
which is boundary ATT with Neumann boundary conditions and function (f)^ {y) = 
J-'ooiy) can be calculated using the perturbation theory. It can be written in the form: 



-E 



-2m,i?/ cosh 9 



AVh cosh2^-cos2(4/2) 



(9^ 



Here //j are the fusion constants and 6'-j are the fusion angles. The constants //j are 
different from zero and equal to 1 or —1 only if particle j is a bound state of two 
particles i. In this case the fusion angle is defined by relation: mj = 2mj cos(6'^j/2). 
The complete list of constants //j and angles ej^ can be found in Ref. |2^. The con- 
stants dj describing one particle contributions are different from zero exactly for the 
same particles j that and the constants dj calculated in previous section. For these 
particles they can be expressed in terms of residues of amplitudes Rj{e) (see Eq.(|65D) 
as: 

TT 



(99) 



As an example of the application of Eq. (PSD we give the explicit expression for 
function (j)^ {y) in D^. ATTs: 



'{y) 



vr 



[r/2]-l 



r de 

e-^-- 1^(91 -9r-i-^ + f^[ 



(100) 



where 



. Nexp(-4v/2sin(7r///i)m?/cosh( 
giie) = sm(27r///i)- 



cosh e — CO?? {t:1 / h) 



cosh e — sm. [2tiIi h) 
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A plot of the functions (j)Q = — cq ■ 4> and (pQ = —Cq ■ (j)^ for ATT is given in Fig.4. 



It is straightforward to derive from equations ( |100| , |101| ) that 



If we take into account Eg. (|77|) we can easily find that boundary value 0^(0) coincides 
with vector 9^ defined by Eq.(^). The same result can be obtained for all simply 
laced ATTs. 

2. In this paper we derived the weak and strong coupling limits of function J-'b{y)- 
In the intermediate region the long distance behavior of this function can be expressed 
in terms of boundary amplitudes Rj and form factors . Some of form factors in 
ATTs are calculated in Ref. ||10||. At short distances we can calculate J-'b{0) which is 



equal to vector 9. Using boundary conditions and equations of motion we can express 
first two derivatives of this function through the boundary one-point function Gb{q) '■ 

= ^^j:e,GB{he,)- ^:(0) = 6VEe,G^(2fee,). (103) 
^ i=o i=o 

To calculate the next derivatives of this function we need the boundary VEVs of the 
descendant fields. At present the problem of their calculation is not solved yet. We 
hope to return to this problem in the next publications. 

3. In the main part of this paper we studied boundary simply laced ATTs with 
dual to Neumann boundary conditions (p^). At the end of this section we consider the 
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classical limit of non-simply laced ATTs with the same boundary conditions. These 
theories are dual to the ATTs with dual affine algebras and boundary conditions 
(PI). In the non-simply laced case the corresponding boundary vacuum solutions 
(f){y) satisfy the same Eqs.(|35D with modified boundary conditions at y = which can 
be written as: 



dy(f) = mJ2y'^ni/(^hiew{ei- (f))- (104) 

The solutions (pnsi of these equations can be obtained from the solutions (psi for 
simply laced ATTs described in the previous section by the reduction with respect 
to the symmetries of extended DE Dynkin diagrams (which are the symmetries of 
solutions (psi)- This reduction can be described in the following way. Afiine Dynkin 
diagrams for non-simply laced algebras can be obtained by folding the nodes of di- 
agrams for simply laced algebras G. Let us the set of nodes {j} of affine Dynkin 
diagram for G corresponds after the folding to the node j of the diagram for non- 
simply laced algebra Q. We denote as ej the root of G with arbitrary J G {j} and as 
ej the corresponding root of Q. Then it holds: 

ej ■ (t)nsi{y) = ej ■ (psiiv) (105) 

This relation defines completely the solutions of Eqs.( |55| , [104D for all for non-simply 
laced algebras Q. The ratio of corresponding boundary ground state energies is equal 
to the ratio of Coxeter numbers of Q and G. 

Extended Dynkin diagram for the Lie algebra Br can be obtained by folding nodes 
r and r -|- 1 of the diagram for -Dr+i- In this case 8^'^^\Br) = £^'^''\Dr+i). In the dual 
limit Br ATT with boundary conditions ( |104| ) is described by weakly coupled B^ (or 
^2ri-i) ATT with Neumann boundary conditions. The perturbative calculations and 
analysis of boundary S-matrix for this theory gives us the reasons to conjecture that 
in quantum case the boundary ground state energies for Br and -Dr+i ATTs have a 
similar form: 

(J,. mcos(7r/4-7r/2/i) 

^bound[J^r) - ^ gi^(^^/2/l) COs(7r(l - x) /2h) ^ ^ 

where h now depends on x: h = 2r — x. 

Folding the nodes r + 1 with r + 2 and with 1 of -Dr-+2 diagram we obtain the 
diagram for or ( -D,;._^\) ATT. The similar consideration of this theory leads us to 



conjecture that Sioundi^r) ^I^o described by Eq. ( |106|) where now h = 2r-|-2(l 



■X) 



Folding the nodes k and 2r — k oi affine D2r diagram we obtain the classical 
solution for 5^ boundary ATT. In this case £^^^\B'^) = £^^^\D2r)/2. The solution 
for BGr (or Ag) ATT can be obtained from D2r+2 solution by folding D2r+2 diagram 
with respect to all symmetries. In this case all r particles contribute to the solution 
and £^^^\BCr) = 8^^^\D 2r +2) 1'^- 

The solution for G2 ATT can be obtained from solution by folding the nodes 
1, 3, 4 of affine Dynkin diagram. In this case only the heavy particle m2 contributes 
to the solution and S^'^^\G2) = S^'^^\Di). Contrary to that only the light particle mi 
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contributes to the solution for i^'^^) ATT. This solution can be obtained by 

the reduction with respect threefold symmetry of Eq extended Dynkin diagram. In 
this case S^'^^^G'^) = S'^^^\Eq) /?>. Folding the nodes 1 with 1 and 3 with 3 of the same 
diagram we obtain "one particle" solution for F4 ATT. The "three particle" solution 
for or {Ef^) ATT can be derived by the reduction with respect to Z2 symmetry 
of E^ extended Dynkin diagram. As a result £^^'-\F^) = £'-''^\Er)/2. The solution 
for the last ATT corresponding to the Lie algebra Cr vanishes: (f){y) = 0. 

In this paper we found explicitly all solutions of Eqs.(P^,|l(J4|). However some 



of non-simply laced ATTs require more general integrable boundary conditions [|I2 
which form the parametric families. We intend to describe the corresponding bound- 
ary one-point functions, scattering theory and boundary solutions in a separate pub- 
lication. 
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Appendix A 

In this appendix we give the explicit form for the function T2 in Ei ATT. This 
function can be obtained from the function tq (see Eqs.( P(]|pTD using Eq. ( P^D and is 



(2) (2) 
determined by the following coefficients Y^^^^, a,b,c = 0, 1, 2; "Kqqq 



^2 = T.T.T.yl^cit2Z2nUZ,)\trZrr (107) 



a=0 b=0 c=0 

where 

= = 4cos(27r/9); Fqi'd = = -4cos(7r/9); 
Yool = = -4 cos(47r/9); f/S = V cos2(7r/9); = "V cos2(27r/9); 
_ _4cos2(47r/9)/cos2(7r/9); Y^^} = 128 cos^(47r/9)/ cos2(7r/9); 

Y200 = yh >o?o = yh ySI = yh y^II = yi2y?yi/y2; 
Y120 = z/i2i/f ^2/^/1; Y^qI = yisys^yi/ys; Y^^i = z/isi/f Vs/yi; 
^021 = i/23i/f ^1/2/1/3; >"o?2 = y23y?y3/y2; Y^^o = 2/12; ^202 = 2/13; 

^022 = 2/23; "^2?! = 2/i22/i3>"oii/ (2/22/3); ySI = 2/122/23 >^ioi/ (2/i2/3); 

Yiil = 2/i32/23yi?o/(2/i2/2); Y§i = 2/i232/i22/fV2/3; >^2?2 = 2/i232/i32/fV2/2; 

^122 = 2/1232/232/f V2/1; >222 = 2/?23- 
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Appendix B 

In this appendix we list the coefficients Y^'^l^^^^^ in the Eq.(^) for tq function in 
ATT. It is convenient to simplify the notations. The variables CTj take the values 
or 1. We can introduce the variable 1 < ^ < 4 indicating o"j which are equal to 
one. The corresponding coefficients We denote as Uij^ Vijk 

and ?/i234- For 

example, ?/2 = ^oito; Vis = ^io w and s o on. Then YoSjo = 1 and 

„, _ 43 19 , / 4289 9362" . „, _ 44 , 83 a [298 T58". 

- - - + V ~ ~ Vs' - -- + 75 - 4 ^— - 

„ _ 71 101 I A TSOT 14342 . „ _ 25 19 , c fTj 38~. 

- Y - 7j + y — - ITT' ^23 - - - 7g + ^y - - iTs' 

^, _ 52 I 43 /I /T25 278 . ^, _ 43 , 19 , / 4289 9362" . 

2/24 - -y + 75 - 4y — - 2/34 - -- + + y — - 

2/123 = 2/l2?/l3?/23/(?/l?/2?/3); 2/124 = 2/l22/l4?/24/ (?/l2/22/4) ; 
2/134 = 2/132/142/34/(2/12/32/4); 2/234 = 2/232/242/34/(2/22/32/4); 
2/1234 = 2/122/132/232/142/242/34/(2/12/22/32/4)^- 



Appendix C 

The classical vacuum solutions of Sec.^ have been cross-checked numerically. We 
used Mathematica[]'s routine NDSolve which allow good accuracy control through the 
parameters Pre c, AccuracyGoal , WorkingPrecision. The code for the solution of 



boundary problem ( |35| , |36D has been developed for all cases considered in this paper; 
it is available at [http : //www, f is .unipr . it/^onof ri| . 

The inherent instability of the boundary value problem (i.e. the existence of 
solutions with exponential growth at infinity) makes the direct numerical problem 
rather ill-conditioned; however this fact is completely under control. The numerical 
solutions agree with the exact ones until the unstable modes eventually take over 
by amplifying the unavoidable truncation and algorithmic errors. In Fig. 5 we may 
see that the discrepancy between exact and numerical solution eventually grows with 
exactly the slope dictated by the square root of the largest eigenvalue of the mass 
matrix (the pattern is replicated systematically for all the other Toda systems). 

In case an exact solution were not available, a possible strategy, is given by 
"backward-integration" : starting from large values of y and assuming a general linear 
combination (p = J^dj^j exp{—mi'j y} one can integrate back to 2/ = and fit the 
coefficients dj to match the boundary conditions. This technique is not affected by 
the unstable modes and can give rather accurate values for the asymptotic data. 



^Copyright ©Wolfram Research Inc. 
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DgATT 




Figure 5: The exponential growth of the difference between exact and numerical solution: 
the slope is given by the square root of the highest eigenvalue of the mass matrix. 
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